function phi = phi1(spectrum,m,n,x,y)
LatticeSize = spectrum.LatticeSize;
XYmax = spectrum.XYmax;
phi = 2.0*sin(2*m*pi*x/LatticeSize)*sin(2*n*pi*y/LatticeSize)/XYmax;